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FIGURE  I 1-2 


ILLUSTRATION  OF  AREAS  INVOLVED 
IN  DETERMINING  GRADING  CRITERIA 


Discussion  of  a^,  a^,  and 

Throughout  the  presentation  of  this  report,  the  quantities  a , a , and  o 

1 M 

are  used  to  characterize  the  ground  surface.  The  units  of  these  variances 
are  wavelengths,  with  one  wavelength  (A)  for  glide  slope  analysis  correspond- 
ing to  approximately  3 feet.  Therefore,  a value  of  o * 0.  5 = 0.  5X  translates 
as  a = 1.5  feet. 

To  properly  understand  the  concept  of  variance,  it  is  useful  to  consider 
how  one  might  characterize  the  ground  surface.  It  is  customary  to  define  the 
mean  height.  Ah,  as  zero  with  the  ground,  then  varying  plus  or  minus  about 
this  zero  level.  Clearly  the  only  remaining  quantity  to  describe  the  surface 
is  the  spread  in  the  distribution  of  height.  Consider  figure  II-Al,  where  the 
distribution  of  surface  heights.  Ah,  is  illustrated.  The  area  under  the  curve 
between  ±a  corresponds  to  about  70  percent  of  the  total  ground  surface.  This 
means  that  when  one  specifies  a o = 0.  5 = 0.  5X  or  o = 1.5  feet  for  glide  slope 
systems,  70  percent  of  the  ground  lies  between  ±1.  5 feet  about  some  zero 
level  drawn  through  the  surface. 


Figure  II-Al  Distribution  of  Heights,  Ah,  and  Illustration  of  a 
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1.  INTRODUCTION 


Scattering  of  electromagnetic  waves  from  rough  surfaces  has  been  studied 

through  the  Kirchoff  approach^  for  very  rough  surfaces,  the  small  perturba- 

2 3 4 

tion  method  of  Rice  for  slightly  rough  surfaces,  and  the  composite  model  * 

for  surfaces  with  both  small  and  large  undulations.  * Surfaces  characterized 

by  more  than  c.ie  statistically  distributed  roughness  have  been  discussed  by 

6 7 

Beckmann.  The  approach  used  in  the  present  work  is  discussed  by  Weiss 

and  consists  of  treating  the  specular  and  diffuse  components  of  the  scattered 

field  separately.  The  diffuse  scattered  rsidiation  is  modeled  on  the  basis  of 

subdividing  the  relevant  scattering  area  into  smadl  rectangular  blocks  and 

using  the  plane  wave  approximation  to  get  the  scattered  fields.  For  the 

g 

specular  portion  of  the  radiation,  the  half-plane  solutions  of  Senior,  Woods, 
and  Bromwich  are  modified  to  yield  expressions  for  scattering  from 
slightly  rough  half-planes.  The  results  are  very  similar  to  the  conventional 
Rayleigh  result  for  specular  scatter  (reference  7,  p.  5). 

Once  expressions  for  the  scattering  of  specular  and  diffuse  radiations  are 
obtained,  they  are  applied  to  the  problem  of  determining  aillowable  grading 
for  instrument  landing  systems.  Current  grading  of  ILS  sites  requires  rathe 
stringent  preparation  of  the  ground,  and  the  application  of  this  paper  is  to 
develop  grading  criteria  that  reflect  the  statistical  nature  of  the  disturbance 
in  observed  signal  rather  than  simply  to  consider  the  maximum  allowable 
phase  differences  as  the  basis  for  grading  criteria,  as  is  now  done. 

The  results  of  this  work  as  will  be  seen,  indicate  that  the  current  grading 
is  too  stringent  and  that  substantial  cost  savings  can  be  effected  through  a 


♦Numbered  references  appear  at  the  end  of  this  document. 


more  systematic  approach.  In  specifying  the  grading  criteria,  the  expressions 
developed  allow  for  the  treatment  of  several  statistically  different  surfaces 
with  the  requirement  that  the  solutions  be  self-consistent  when  the  roughness 
of  each  adjacent  surface  is  equal. 


2.  SPECULAR  CONTRIBUTION  OF  SCATTERED  FIELD 


In  developing  expressions  for  the  field  scattered  from  a slightly  rough 
half-plane,  it  is  appropriate  to  consider  the  two-dimensional  dual-integral 
formulation  of  the  problem^  ^ and  then  generalize  to  the  three-dimensional 
case  as  outlined  by  Senior.  Consider  the  plane  electromagnetic  wave  with 
harmonic  time  dependence  exp(i<A>t)  incident  obliquely  at  an  angle  on  a 
slightly  rough  half- plane,  z = Ah(x,  y),  and  x < x^,  where  <Ah>=  0 as  illustra- 
ted in  figure  II- 1.  If  XYZ  defines  a set  of  space  axes  with  unit  vectors  u J» 
and  ic,  a local  coordinate  system  on  the  surface  can.  be  defined  with  axes 
X'Y'Z'and  unit  vectors  andic,'  respectively.  For  slightly  rough  ground, 

the  variation  of  Ah(x,  y)  with  x and  y is  small  compared  to  unity,  and  if  the 
variation  in  y is  small  compared  to  x,  the  surface  unit  vectors  are  related  to 
the  space  unit  vectors  through 

‘i'  = ^ ^ 

(lb) 


Ah  + 

X 


where  Ah^  = 3Ah/3x.  The  incident  plane  wave  satisfies 

» 

-ilw  -ik  (x  cos  - z sin  a.)  (2) 

e = e 0 0 

-ik  [x(cos  a'  + Ah  sin  a')  - (z  - Ah) (sin  a'  - Ah  cos  a')l 
= e 0x0  0x0 

where  is  the  incident  angle  defined  relative  to  the  surface  local  coordinate 
system.  The  dual  integral  equations  describing  the  scattered  fields  at  a point 
on  the  local  surface  (z  = Ah)  because 


1 - 


ik  (x  -x')  (M-'-Ah  >yi-P'  ) , , ik  (x  -x')  (p'  + Ah  Vl-p'  ) {3a) 
e e X dp'  = -e  e 0 x 


rcsaasss 


ik(x  -x')(Ki'-Ah  Vl-fA'”)  , , » , 1 

P(HL')e  e X ' d»A'  = 0 (x>x  ) 

e 


(3b) 


where  K'  = cos  = cos  and  P(ja')  is  an  appropriate  Fourier  weighting 

coefficient  defined  by  (see  reference  1,  p.  567) 


P(p-)  = ^ 


(4) 


Equations  (3a)  and  (3b)  can  be  further  simplified  if  the  variations  across  the 
surface  are  characterized  by  large  correlation  distances  (small ^h^).  Here 
the  slope  terms  can  be  ignored  and  the  scattered  field  simply  becomes  that 
for  the  usual  half-plane  solutions  with  z replaced  by  z - Ah  (see  reference  ll, 
p.  567).  Generedizing  the  two-dimensional  solution  to  three  dimensions  as 
outlined  by  Senior  (reference  8),  the  scattered  Y-component  of  the  field  from 
a dipole  along  the  Y-axis  becomes 


:^^(Ah)  = - isi.  fdafd^fdy  .8.in(l/2)«  sin(l/2) 
YY  2n  J J J COSa  + cosy 


-ik|(x  -X  )cosa  cosp  + (z  - Ah)sina  cosp  - y sinp) 
* 6 GO  O O 

^ ^ik  l(x^-x)cos  y cosp  + ( z _ Ah)  siny  cosp  +y  sinpi 


(5) 


where  the  integrals  over®,  p,  and  7 are  the  steepest  descent  paths  of  refer- 
ence 2;  x^,  y^,  denotes  the  position  of  the  dipole;  x,  y,  z denotes  the 

observer  point;  and  x denotes  the  position  of  the  half-plane  edge.  The  nota- 

® g 

tion  for  the  field  in  equation  (5)  corresponds  to  a scattered  field  C..»  where 

i denotes  the  field  component  in  space  and  j denotes  the  axis  along  which  the 
1 2 

dipole  is  parallel.  The  dependence  of  the  specular  scattered  field  integrals 

on  the  reflection  coefficient  for  horizontally  polarized  radiation,  R~  (a),  has 

o 

been  avoided  since,  to  a good  approximation  for  horizontally  polarized- rad ia- 

13 

tion  scattered  off  a plane  earth. 


ir;(»)I  =1 


(0  < a < 4 ) 


(6) 


.1 


where  the  bulk  of  the  contribution  to  the  integr^ds  comes  from  the  specular 

angle  with  o < 4°  (^  = 3 m).  For  glide  slope  arrays,  the  dipole  is  parallel  to 

the  Y- space  axis  with  the  aircraft  antenna  similarly  polarized  so  that  the 
S 

field  component,  Cyy’  case  of  interest. 

Expanding  the  exponentials  in  Ah  and  reducing  the  integrals  of  equation  (5) 
to  Fresnel  integrcds  in  the  manner  outlined  by  Senior, 


mn 


m+n 


m!n! 


'0Z  ’ 


m 


(7) 


az'  YY 


where 


I?Y  = I^rIO 


(8) 


and  where  k = Ztt/A 


i;r/4 


(9a) 

(9b) 


where  u = R or  S and 


/krr 

• = zJ-  ° 

u IfR,  + u 


cos(^  A $ )/2 


u = [(x-x  )^  + (y-y  )^  + (z  ± z )^] 
L O O O J 


2II/2 


(10a) 


(10b) 


In  equations  (10a)  and  (10b)  the  plus  sign  corresponds  to  u = S and  the  minus 

sign  to  u = R.  The  angles  6 and  0 are  defined  by  z = r sin  0,  x-x  = r cos  0 

o e 

(x  > X ),  z = r sin  0 , and  x - x = r cos  0 (x  > x ).  The  remaining 
eooo  eoo  oeo  ® 


variables  are 
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The  solutions  to  the  specular  scattered  field  presented  in  equation  (7)  is 
equivalent  to  the  Rayleigh  (see  reference  5,  p.  240)  result  for  rough  scatter 
except  that  it  has  been  developed  for  a slightly  rough  half- plane.  This  solu- 
tion is  quite  adequate  for  treatment  of  glide  slope  problems  where  the  specular 

o 14 

angles  are  always  small  (<4  ).  It  was  found  that,  for  ground  roughness  with 
variances  a < 1.  7^  (assuming  normally  distributed  ground),  a second-order 
expansion  in  equation  (7)  is  sufficient. 
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3.  THE  DIFFUSE  CONTRIBUTION  TO  THE  SCATTERED  FIELDS 

The  diffuse  scattering  contribution  to  the  fields  is  not  easily  determined 
from  theoretical  considerations.  In  general  this  contribution  is  linearly  in- 
dependent of  the  specular  scattered  fields  and  can  be  written  in  the  form  (for 
the  Y- component  from  a dipole  along  the  Y-axis) 


■ DIF 
YY 


(12) 


where  « = a (a)  is  the  incident  angle  measured  from  the  surface  to  the  incident 

wave  vector,  R is  the  reflection  coefficient,  and  p_  is  the  diffuse  scattering 
o D ° 

coefficient.  The  reflection  coefficient,  R^,  is  that  for  horizontally  polar- 
ized radiation,  and  since  the  integral  in  equation  (12)  includes  angles  where 
R^  is  appreciably  different  from  unity,  it  must  be  included  in  the  expression 

for  diffuse  scatter.  The  functional  behavior  of  R (a)  was  obtained  from 

o 

Beckman  and  Spizzichino  (reference  13,  p.  220)  as  was  the  behavior  of  Pq(“) 
(reference  13,  p.  340).  The  diffuse  scattering  coefficient  functionally  depends 
on  the  parameter 

4<r.^  sin  a 


« = 


(13) 


where  is  the  variance  of  the  roughness. 

The  integrcd  appearing  in  equation  (12)  was  evaluated  numerically  assuming 
the  surface  to  be  subdivided  into  smedl  rectangles  subject  to  the  condition 
that  the  incident  radiation  could  be  approximated  as  a plane  wave  over  each 
small  rectangle.  This  criterion  was  used  to  determine  the  appropriate  size 
of  each  rectangular  subsection  of  the  surface.  The  scattered  field  due  to  each 
rectangular  section  was  calculated  within  the  framework  of  the  vector  Kirchoff 
integral  representation  of  diffracted  radiation,  with  the  result 
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EYY^(a)=  - 


ik 


ikr 


Ztt 


s A 


ik(r  + r.) 
e » A »i>2 


sin  A sin  B 


(14) 


B 


where  ^ is  the  direction  of  the  dipole  and  measured  scattered  polarization, 

is  ilie  “nit  scattered  wave  vector,  is  a unit  vector  normal  to  the  surface 

element,  R is  as  defined  by  equation  (10b),  r is  the  distance  from  the  plate 

s 

center  to  the  source,  r is  the  disteUice  from  the  plate  center  to  the  observer, 
a^  and  a^  are  the  length  and  width  of  the  rectangular  section,  and 

A = k cos  (/S/2)  sintf  coa<t>  a^ 

B = k cos  (j8/2)  sintf  sin^  a^ 

P = angle  between  ^ and  ft 

where  6 and  4>  are  spherical  coordinates  defining  the  direction  of  a vector  ft^ 
in  the  plane  of  j and  ft^  that  bisects  the  angle  p.  The  limits  on  the  surface 
of  integration,  S,  of  equation  (12)  were  taken  as  that  v«d.ue  of  £ such  that 
was  < 0.  01. 
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If  p defines  the  position  of  the  observer  (p(x,  y,  z))  and  p'(m)  defines  the 
position  of  the  mth  dipole  in  the  radiating  array,  the  difference  in  depth  of 
modvilation  (DDM)  signal  (see  reference  12)  is  defined  by 


DDM(p)  = NR 


^n  jSBO  _ , ,/  u 

Eyy(P.  P'(in)) 
m 

En  jCSB  _ , 1/  w 

1 


SBO 

where  N is  an  appropriate  normalization,  d is  the  dipole  amplitude  and 

CSB 

phase  of  the  mth  dipole  radiating  the  sideband  signal,  and  d^  is  the  dipole 
and  amplitude  and  phase  of  the  mth  dipole  radiating  the  carrier  signal.  The 
complex  field,  Eyy»  contains  both  the  scattered  and  direct  radiation.  The 
normalization,  N,  is  obtained  by  requiring  the  DDM  to  be  ±75  at  angles 
of  ±0.  35°  about  the  desired  aircraft  glidepath. 

Contributions  to  the  DDM  from  the  diffuse  and  specular  radiation  are 
statistically  uncorrelated;  consequently,  it  is  useful  to  treat  each  type  of 
radiation  separately,  with  the  requirement  that  the  total  DDM  (being  the  sum 
for  uncorrelated  components)  be  less  than  a prescribed  amount  that  is  based 


on  the  tolerances  allowed  for  DDM  deviations, 
becomes 


The  total  diffuse  field 


„ _ / -ikR 

-k2(e 

^YY  ■ V kR 
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■which  yields  for  the  diffuse  component  of  the  DDM 


D 


DDM(p)  = NR 


in 

n CSB  D 

^ 1 ^YY 

m 


(17) 


The  average  over  roughness  in  the  diffuse  component  is  implicity  contained 
in  the  expression  for  where  the  dependence  is  expressed  in  terms  of 
the  variance  o^^. 

Evaluation  of  the  contribution  to  the  specular  component  of  the  DDM 
requires,  however,  that  the  signal  be  averaged  over  the  roughness  function 

g 

Ah  = Ah(x,  y).  In  addition  to  finding  the  average  DDM  (specular  contribu- 

s s 

tion),  it  is  necessary  to  examine  the  variance  in  specular  DDM  , o $ 

DDM 

since  this  quantity  is  a measure  of  the  spread  of  the  specular  component 

and,  consequently,  it  represents  the  true  measure  of  the  deviation  about 

the  mean  of  the  specular  component.  The  distinction  between  the  two  pro- 

cedures  for  calctilating  the  components  of  DDM,  DDM  , and  DDM  is 

simplified  if  one  recognizes  that,  while  equation  (17)  represents  a worst 

case  calculation  for  the  diffuse  contribution  owing  to  the  fact  that 

implicitly  contains  the  roughness  already  averaged,  the  expressions  to  be 
g 

developed  for  constitutes  the  worst  case  for  specular  DDM  about 

the  mean,  <DDM®>. 

C 

In  evaluating  the  average  specvilar  DDM,  < DDM  >,  and  the  variance, 

S 16 

^Ddm’  desirable  to  introduce  the  bivariate  distribution. 


f(Ah,,  Ah.)  = 
1 J 


1 


Zrra.a.  Vl-pf. 
I J ij 


■2( 


B:-)  • (■^)  ■ -'-.{-tm 


(18) 
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where  o.  and  O.  are  the  variances  in  height  of  Ah.  = Ah(x.,  y.)  and  Ah.  = 

1 J 1 1 1 j 

Ah(Xj,  yj.  The  correlation  is  measured  by 


■Ah.  - Ah. 
1 1 


P..  = e 


For  the  cases  considered  here,  it  was  found  that  p.^  was  vanishingly 

small.  It  was  found  that  the  dependence  of  <DDM^  > and  on  the 

^ DDM 

correlation  coefficient  was  insignificant,  and  appropriate  expressions  for 
these  quantities  to  the  second  order  are 


<DDM  > = NR 


["o  1 * '=„  J (b  ^ 

o 

7 


2 b .1 
> J 2jJ 


(20a) 


where 


„n  ,SBO  g-i^S(j)  g ^ -| 

" ^i  *^3  * ^ L kR(j)  ■ KS(j)  ^ ^YY 


(ZOb) 


(21a) 


,SBO 
a. . = d. 
y J 


r 2 9 2 9 9 ) S 1 

r lYY(R(j).S(j))J  (21b) 
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Grading  criteria  are  based  on  subdividing  the  ground  in  front  of  the  antenna 
into  the  three  principal  regions  of  figure  II-2.  A further  development  for  the 
specular  DDM  is  thus  required;  namely,  it  must  provide  for  two  areas 
characterized  by  different  roughness,  and  Applying  the  principle 
of  superposition  to  the  half-plane  solution  described  earlier,  it  is  possible 
to  characterize  the  ground  of  figure  II-2  by  a half-plane  of  width  (-<»  to  x. ) 

3 

and  a strip  (Xj^  to  x^).  Any  expression  for  the  resulting  variance  must 

be  self-consistent  in  that  when  = er^  solution  is  the  same  as  for  a 
half-plane  of  width  (-<»  to  x^  + x^).  Such  a development  can  be  obtained  based 
on  the  earlier  discussion  with  the  result 

.22  C,.Cl.«,,c^p2^-a^)-2  2,C,.)^V2'‘l2 

2 2. 


" f Sj '4  - 4']  - «’,.<'2»' 


(23) 
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I t I 

Here  the  quantities  and  are  defined  in  the  saxne  way  as  for 

the  C s appearing  in  equations  (24),  except  that  x = x . The  correlation 
coefficient  now  becomes 


Pl2  = ® (25) 

It  is  clear  that,  when  = O^,  equation  (23)  is  self-consistent.  An  expression 
for  <DDM^(a^,  0^)^  is  easily  obtained  from  equation  (20a). 


IS 


5.  CALCULATION  OF  CORRELATION  DISTANCES 


In  figure  2 the  areas  characterized  by  and  also  included  a dependence 
on  the  parameters  1 and  1 . The  quantities,  the  correlation  distances, 

X Ct 

specify  the  distance  between  peaks  in  the  roughness  and  in  general  depend 
only  on  the  statistical  value  of  the  ground  (here  assumed  random).  To 
simplify  the  calculation  of  roughness,  it  has  been  assximed  that  the  slope 
terms  appearing  in  equations  (3a)  and  (3b)  were  ignored.  Clearly,  deter- 
mination of  an  appropriate  correlation  distance  should  be  consistent  with 

this  approximation. 

7 

Bareck  has  found  that  the  density  of  specvilar  points  satisfies 


7.  255 

”a  ~ 2,2 

TT  Jt 


tan^> 


whe  re 


tan7  = 


2 2 
sin  9.-2  sin9.  sin  9 cos  9 + sin  9 
i 1 s s s 

cos  9,  + cos  9 
i s 


2 2 2 

with  i being  the  correlation  distance,  s = 4o  /i  is  the  mean  slope  at  a 

point  on  the  surface,  9.  is  the  incident  angle,  and  9 is  the  scattered  angle. 

1 s 

The  total  number  of  specvilar  points  contained  within  an  area  satisfies 


Integration  of  equation  (28)  numerically  yields  the  results  for  the  first  area 

(characterized  by  a,)  that  N.  = 1 when 

1 A 
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i(X) 

!l 

n 

0.  1 

2 

0. 

05 

1 

* 

0.  2 

3.5 

0. 

057 

0.3 

5 

0. 

06 

• 

0.4 

6.5 

0. 

06 

Since  for  the  sirall  specular  grazing  angles  appropriate  to  glide  slope 
systems  = 1,  it  is  clear  that  slopes  of  about  0,  05  are  required  to 
satisfy  the  approximations  used  in  this  calctilation.  The  values  of  i(X) 
above  indicate  how  the  ground  correlations  must  be  graded  for  various 
choices  of  o^.  A similar  correlation  for  with  0,  01  yields  slopes 
of  0.  02.  These  criteria  will  be  assumed  when  calculating  the  requirements 
on  grading. 


i 
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6.  NUMERICAL  RESULTS  FOR  GRADING  CRITERIA 


Figures  II-3a,  b,  c,  d,  e Illustrate  the  effect  of  the  diffuse  scattered  field 
on  the  DDM  for  a standard  earth  and  a surface  infinite  conductivity.  These 
figures  are  for  capture  effect  arrays  and  represent  the  error  in  DDM,  DDM^, 
compared  with  infinite  ground  results  for  range  points  at  30000,  15000,  10000, 
5000,  and  3000  feet  from  threshold.  These  points  were  chosen  as  repre- 
sentative of  the  flight  path.  The  array  is  located  1100  feet  behind  threshold 
and  300  feet  left  of  runway  centerline.  The  worst  derogation  is  seen  to 
occur  at  30000  feet  from  threshold.  Table  11-2  illustrates  the  relative  values 
of  versus  the  error  introduced  into  the  DDM.  It  is  clear,  for  example, 

that  values  of  = 0.  032X  for  metal  surfaces  and  9.  = 0.  044X  for  standard 
D D 

earth  correspond  to  errors  on  the  order  of  10  pA  in  the  observed  DDM  for 
capture  effect  systems. 

Figures  II-4a,  b,  c,  d,  e and  II-5a,  b,  c,  d,  e show  similar  plots  for  the 
nvill  reference  and  sideband  reference  systems.  Again  the  cases  of  a metal 
surface  and  standard  earth  are  plotted.  Table  II-2  illustrates  the  range  of 
values  of  for  these  systems  versus  maximum  observed  error  in  DDlvt,  It 
is  clear  that  the  null  reference  system  performs  best  with  regard  to  ground 
roughness  and  diffuse  scatter  followed  by  the  capture  effect  and  sideband 
reference  systems.  This  result  is  in  keeping  with  the  ground  current 
structure  near  each  array:  The  sideband  reference  has  the  highest,  followed 
by  the  capture  effect,  and  the  null  reference  has  the  lowest  ground  currents 
near  the  antenna.  Clearly,  the  higher  the  ground  currents  in  the  diffuse 
region,  the  greater  the  possibility  for  signal  derogation  due  to  diffuse 
scatter  from  rough  ground.  In  addition,  the  standard  earth  case  demonstrates 
substantially  reduced  error  over  a flat  metal  surface.  The  resiilts 
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r 


for  the  metal  surface  are  qualitatively  in  agreement  with  the  experimental 

1 8 

observations  of  Lucas.  Conventional  grading  of  ILS  sites  is  to  within 

15 

±0.  1 foot  during  the  first  1000  feet  of  area.  This  wl  aid  correspond  to 
a surface  roughness  of  about  ±0.  03A.,  which  indict*  ss  that  the  observed 
error  due  to  diffuse  scatter  would  be  less  than  2.  5 |jlA,  7.  5 pA,  and  10  pA 
for  null  reference,  capture  effect,  and  sideband  reference  systems, 
respectively. 

Figure  II-2  shows  a semicircular  area  of  radius  r^  representing  the 

region  of  dominant  contribution  to  diffuse  scatter.  It  is  clear  that  r^  is 

determined  by  the  point  at  which  R (a)  reach  some  minimum  value  so 

o D 

that  the  remaining  contribution  to  diffuse  scatter  is  negligible.  Defining 
this  point  as  the  one  at  which  95  percent  of  the  diffuse  radiation  is  included 
in  the  expressions  for  diffuse  field,  one  can  specify  a grading  limit  on  the 
radius  of  this  diffuse  region  as  a function  of  the  roughness  o^.  Table  II-3 
illustrates  the  results  for  this  radius  versus  a^.  As  increases,  more 
diffuse  scatter  is  observed  coming  from  segments  of  the  ground  located 
further  from  the  antenna,  and,  consequently,  the  amount  of  ground  requiring 
grading  to  minimize  diffuse  scatter  increases.  It  is  important  to  note 
that,  although  the  sideband  reference  has  the  highest  ground  currents  near 
the  antenna,  they  occur  at  closer  distances  to  the  array  since  the  antenna 
elements  are  closer  to  the  ground  surface;  consequently,  a smaller  segment 
of  ground  under  this  array  is  subject  to  diffuse  scatter. 

Figure  II-6a,  b,  c,  d,  e,  figure  II-7a,  b,  c,  d,  e,  and  figure  II-8a,  b,  c,  d,  e 
illustrate  the  specular  contribution  to  the  variance  in  the  DDM  versus 
and  (the  two  roughness  factors  characterizing  the  areas  in  front  of 
the  arrays  as  shown  in  figure  II-2).  Several  effects  are  clear  from  these 
figures.  The  contribution  from  the  first  area,  o^,  dominates  the  behavior 
of  the  variance  in  DDM,  and  only  at  small  and  large  distances  from  the 
array  does  the  second  area,  characterized  by  introduce  a significant 
effect.  This  is  because,  at  small  distances  from  the  array,  the  edge  of 
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the  graded  area  that  wae  taken  to  be  2000  feet  noticeably  degrades  the 
pattern.  This  effect  Is  Illustrated  In  figures  II-9,  II-IO,  and  II-ll,  where 
the  Influence  of  a truncated  half-plane  edge  on  the  DDM  Is  plotted  for  each 
system.  At  very  large  distances  from  threshold,  substantial  specular 
scatter  occurs  from  the  second  region  owing  to  the  grazing  angle  of  the 
reflected  radiation,  and  one  again  expects  o^  to  become  important. 

Comparison  of  the  behavior  of  the  three  systems  demonstrates  that  the 
sideband  reference  performs  most  poorly,  followed  by  the  capture  effect, 
with  that  followed  closely  by  the  null  reference.  This  is  in  keeping  with 
the  fact  that  the  dominant  roughness  effects  come  from  the  area  closest 
to  the  array  and  that  the  ground  currents  are  highest  for  the  sideband 
reference,  followed  by  the  capture  effect  and  null  reference  arrays  in 
this  area.  It  should  be  realized  that  the  variance  in  the  DDM  is  being 
plotted,  and  in  effect  this  represents  the  allowable  spread  in  DDM  that 
could  occur.  (It  accounts  for  the  spread  99  percent  of  the  time.  ) 

Calculation  of  the  probable  limits  of  signal  derogation  must  include  the 
variance  due  to  specular  scatter,  the  change  in  the  signal  position  about 
the  infinite  plane  case  (figures  11-9,11-10,  and  II-ll),  and  the  contribution 
from  the  diffuse  scatter.  The  area  to  be  graded  is  illustrated  in  figure  II-2, 
where  the  x-distances  were  determined  by  existing  criteria  and  the  y- 
distances  were  determined  by  the  point  at  which  lateral  edges  yield  no 
substantial  effect  (see  reference  12).  The  diffuse  region  is  specified 
as  in  table  II-2. 


7.  DISCUSSION 


A technique  for  describing  specular  and  diffuse  scatter  from  slightly 
rough  half-planes  has  been  developed.  The  resvilts  are  shown  to  be  applicable 
to  areas  characterized  by  varying  ground  roughness  statistics.  Application 
of  that  technique  for  describing  ground  scatter  to  ILS  sites  demonstrates 
that  the  existing  requirements  on  overall  grading  are  presently  too  stringent. 

The  results  for  diffuse  scatter  indicate  that  current  grading  is  sufficient 
to  minimize  the  effects  of  such  scatter  and  to  yield  patterns  within  current 
ILS  glide  slope  requirements.  The  diffuse  scatter  resvilts  are  in  qualitative 
agreement  with  existing  model  range  results  and  practical  observations.  The 
derogation  due  to  specular  scatter  indicates  that  a substantial  relaxation  in 
grading  is  possible.  Figure  11-12  represents  a typical  cost-effectiveness 
plot  versus  where  the  allowable  diffuse  derogation  of  2.  5 ^A,  7.  5 ixA, 
and  10  nA  for  null  reference,  capture  effect,  and  sideband  reference  has 
been  used.  As  the  figure  indicates,  the  proposed  cost  of  site  grading  is 
considerably  less  than  current  costs. 
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